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USAFSAM  ECG/VCG  DIGITIZING  AND  AVERAGING  SYSTEM 


INTRODUCTION 

Computers  are  being  used  more  and  more  to  aid  physicians  in  analyzing 
electrocardiograms.  This  aid  ranges  from  filtering  the  cardiograms  and 
amplifying  them  to  enhance  their  presentation  to  the  cardiologist  all 
the  way  to  analyzing  the  cardiograms  and  making  diagnoses.  The  first 
step  in  using  the  computer  to  aid  the  cardiologist  entails  converting 
the  cardiogram  analog  signals  into  digital  form.  Next,  prior  to  any 
automated  analyses,  the  digitized  signals  must  be  filtered,  to  remove 
noise,  and  stored  in  the  correct  format. 

At  the  USAF  School  of  Aerospace  Medicine  (USAFSAM) , it  is  common 
practice  in  the  Clinical  Sciences  Division  to  adjust  the  gain  of  the 
instrumentation  amplifiers  and  apply  a d.c.  bias  to  maximize  the  range 
of  the  signal  recorded  on  paper.  This  allows  a larger  picture  for  the 
cardiologist  to  examine.  In  the  Calibration  section,  we  discuss  the 
equipment  used  in  acquiring  data,  the  equipment's  characteristics,  and 
the  procedure  used  for  determining  the  actual  voltage  level  of  the 
digitized  signal.  The  procedure  is  based  on  passing  an  a.c.  reference 
signal  of  known  amplitude  through  the  system  after  the  amplifier  gains 
are  adjusted. 

The  procedures  used  for  removing  noises  from  the  digitized  cardio- 
grams are  discussed  in  the  Digital  Filter  section.  First,  straightening 
the  baseline  and  shifting  to  zero  removes  low-frequency  noise.  After- 
ward, the  noise  which  is  uncorrelated  with  itself  from  beat-to-beat  is 
removed  by  averaging  successive  beats.  The  Digital  Filter  section  also 
contains  the  results  of  verification  tests  that  were  run  to  prove  that 
the  algorithms  were  performing  adequately. 

Later  sections  provide  the  formats  used  in  outputting  the  data,  a 
detailed  set  of  operating  instructions,  and  the  conclusions  derived  from 
this  work  for  future  work. 


CALIBRATION 

ECG/VCG  data  to  be  digitized  and  later  processed  are  supplied  to 
the  signal  conditioning  and  digitizing  system  from  an  analog  instru- 
mentation magnetic-FM-tape  recorder/reproducer.  Prior  to  1976,  the 
signal-to-noise  ratios  from  ECG  signals  reproduced  from  tape  were 
limited  to  40  dB.  New  equipment  being  installed  should  boost  signal-to- 
noise  ratios  to  above  50  dB. 
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Analog  Filter 

Each  of  three  channels  of  ECG  data  (assumed  to  be  three  orthogonal 
leads)  from  the  tape  reproducer  is  first  passed  through  an  analog  4 th- 
order  Butterworth  low-pass  filter,  which  has  a d.c.  gain  of  zero  and  a 
100-Hz  cutoff  frequency.  Within  the  passband  the  filters  are  maximally 
flat,  while  outside  the  passband  their  attenuation  increases  at  a rate 
of  80  dB  per  decade.  These  filter  characteristics  were  chosen  to 
improve  the  signal- to-noise  ratio  from  the  tape  reproducer  and  guarantee 
satisfaction  of  the  Nyquist  sampling  criterion  for  sampling  rates  as  low 
as  200  Hz.  While  tapes  recorded  from  our  VCG  sources  (Hewlett-Packard, 
VCG  Programmer,  model  1507C)  may  contain  information  in  the  0-  to 
200-Hz  band  because  of  the  cutoff  frequencies  of  the  preamplifiers, 
tapes  recorded  from  the  treadmill  ECG  source  (Hewlett-Packard,  Pre- 
amplifier, model  8811)  contain  information  in  only  the  0-  to  100-Hz 
band.  In  addition,  it  is  doubtful  that  surface  ECG  data  contain  sig- 
nificant information  in  frequency  bands  above  100  Hz  [5]. 


Amplifier 

Each  filtered  output  is  amplified  by  a factor  of  3 and  presented  to 
the  input  of  the  analog  multiplexer.  A gain  of  3 was  chosen  to  maximize 
the  quantization  accuracy  without  driving  the  multiplexer  into  its 
nonlinear  region.  The  range  of  the  output  of  the  analog  tape  recorder 
is  ±1.4  V,  while  the  range  of  the  input  to  the  multiplexer  is  ±5  V. 

Multiplexer  and  Digitizer  ' 

The  multiplexer,  model  MP4110,  and  the  analog-to-digital  converter, 
model  MP2812 , are  both  products  of  Analogic  Corporation.  The  converter 
provides  12  bits  of  digital  output  within  an  absolute  accuracy  of  ±0.02% 
of  full  scale.  Total  conversion  time  for  one  channel  is  4 usee.  Word 
output  format  is  2 's-complement  to  match  the  internal  arithmetic  con- 
ventions of  the  Nova  minicomputer/controller. 

Data-Acquisition  Software 

Data  acquisition  was  performed  with  a FORTRAN  callable  assembly 
language  subroutine  which  utilizes  the  Data  General  Real-Time  Disk 
Operating  System  (RDOS) . The  data  were  digitized  and  stored  temporarily 
on  a mass-storage  device.  The  sampling  rate  and  duration,  the  number  of 
channels  to  digitize,  and  the  identification  of  the  output  file  are  all 
defined  in  the  calling  sequence  of  this  subroutine.  These  features 
allow  the  use  of  different  sampling  rates  and  acquisition-period  lengths 
for  the  calibration  and  vectorcardiogram  data.  The  calibration  signal 
was  digitized  at  125  samples/sec  for  6 seconds,  while  the  vectorcardiogram 
was  digitized  at  500  samples/sec  for  12  seconds. 
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Another  feature  of  this  subroutine  is  limit  checking.  To  detect 
overscaling  in  either  direction,  a window  slightly  smaller  than  the 
range  of  the  analog-to-digital  converters  was  defined.  Each  digitized 
value  was  checked.  When  a value  was  outside  this  window,  an  element  of 
an  array  defined  in  the  calling  sequence  for  this  purpose  was  incre- 
mented. This  array  is  assumed  to  two  elements  fo*.  each  channel  digi- 
tized. The  first  of  the  pair  counted  the  number  of  values  exceeding  the 
positive  limit;  and  the  second,  the  number  of  values  below  the  negative 
limit. 


Calibration-Factor  Determination 


After  the  biases  and  amplifier  gains  are  set,  a pulse  train  is 
passed  through  the  system  and  recorded.  The  peak-to-peak  variations  of 
the  pulses  are  0.5  mV,  and  the  pulse  rate  is  30  pulses/min.  These 
pulses  are  used  to  calibrate  the  signals. 

The  calibration  process  requires  two  passes  over  the  data.  The 
first  pass  defines  the  midrange  for  each  channel  as  the  mean  of  the 
maximum  and  minimum  values  for  the  channel.  The  lower  value  of  the 
calibration  waveform  is  determined  by  the  second  pass  and  is  defined  as 
the  mean  of  all  values  below  midrange.  Because  the  analog  filters  cause 
overshoot  on  the  leading  edge  of  the  calibration  pulse  and  because  of 
the  low  duty  cycle,  this  technique  would  not  produce  a reliable  value 
for  the  upper  limit.  To  determine  the  upper  value,  the  distance  between 
midrange  and  maximum  is  divided  into  50  segments  and  the  number  of 
points  in  each  segment  is  determined.  These  frequencies  are  examined, 
and  the  upper  limit  is  defined  as  the  midvalue  of  the  three  consecutive 
segments  whose  sum  frequency  is  the  largest  for  each  channel.  The 
calibration  factor  for  each  channel  is  then  defined  as  twice  the  dif- 
ference of  the  upper  and  lower  limits  for  that  channel. 

DIGITAL  FILTER 

Jointly,  USAFSAM  and  the  Charles  Stark  Draper  Laboratory  have 
developed  a system  which  filters  and  averages  electrocardiograms  and 
vectorcardiograms  (ECG/VCG) . The  system  has  been  coded  in  ANSII  Fortran 
and  implemented  on  a Data  General  Nova  1230  digital  computer.  It 
accepts  three  separate  signals,  each  sampled  at  a rate  of  500  samples/ 
sec,  straightens  and  shifts  the  baseline  to  zero,  and  averages  suc- 
cessive beats  to  produce  a representative  beat  for  each  signal.  The 
slope  of  the  baseline  is  continuously  monitored,  and  beats  corresponding 
to  a baseline  slope  greater  than  2 mV/sec  are  not  used  in  determining 
the  average  beat. 


Slope  Estimation 

At  several  stages  in  the  digital  filter,  curves  are  approximated  by 
a straight  line.  For  instance,  in  the  Baseline  Removal  section,  a 
straight  line  will  be  fitted  to  11  selected  points,  each  80  msec  apart. 
The  midpoint  of  the  straight  line  will  be  used  to  estimate  the  baseline 
of  the  incoming  signal.  In  the  Fiducial  Detection  and  Averaging  section, 
a straight  line  will  be  fitted  to  7 consecutive  samples,  each  2 msec 
apart.  The  slope  of  the  straight  line  will  be  used  as  an  estimate  of 
the  derivative  of  the  incoming  signal  at  the  time  corresponding  to  the 
midpoint  of  the  7 samples.  Finally,  in  the  Computer  Program  Flowchart 
section,  a straight  line  will  be  fitted  to  11  consecutive  points,  each 
80  msec  apart.  The  slope  of  the  line  will  be  used  to  determine  the 
quality  of  the  incoming  data. 

In  all  of  the  cases  mentioned,  straight  liges  are  being  fitted  to 
more  than  two  points;  therefore,  since  two  points  determine  a straight 
line,  the  fit  cannot  be  perfect. 


The  curves  are  approximated  by  the  straight  line  that  minimizes  the 
sum  of  the  square  of  the  error  in  the  fits  at  each  point  (least-square 
error  [9]).  Since  the  same  approach  will  be  used  several  times,  a brief 
description  of  it  follows. 

The  problem  is  to  find  the  straight  line 

x(k)  = b + ks  (1) 


which  minimizes  the  sum  of  the  square  of  the  error  in  the  fits  at  N points 


x ( k ) 


N-l 
2 ' 


N-l 

2 


+ 1 , • . • , — 1,  Of  1,..., 


N-l 

2 

J = [x (k) -x (k) ] 2 


(2) 


The  reasons  for  the  peculiar  and  clumsy-appearing  indexing  of  time 
are  to  simplify  the  calculations  and  to  make  b the  estimate  of  the 
straight  line  at  the  time-interval  midpoint  (k  = O) . The  solution  to 
this  problem  is  [9| 
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Observe  that  the  slope  of  straight  line  s is  the  weighted  sum  of  the 
differences  in  the  last  and  the  first  points,  the  next-to-last  and  the 
second  points,  etc.  Also,  observe  that  b is  simply  the  average  value  of 
the  points . 

Baseline  Removal 

Precisely  defining  the  baseline  is  a very  difficult  task.  Some  say 
that  it  is  the  voltage  measured  between  either  the  T and  P waves  or,  if 
a U wave  is  detected,  between  the  U and  the  P waves.  Others  define  the 
baseline  to  be  somewhere  on  the  P-R  interval,  which  in  many  cases  is 
sloped.  The  problem  is  further  complicated  by  the  fact  that  the  electric 
potential  measured  at  either  of  the  points  fluctuates  as  a result  of 
both  physiological  events  and  the  electrical  instrumentation.  In 
addition,  while  recording  cardiograms,  we  normally  add  a d.c.  bias  and 
adjust  the  amplifier  gain  so  that  the  cardiogram  fills  up  the  allowed 
space  on  a strip  chart. 

In  making  amplitude  measurements , we  use  the  baseline  as  the 
reference;  therefore,  the  computer  must  locate  the  baseline  to  make 
amplitude  measurements.  The  approach  used  with  the  digital  filter  is  to 
define  the  baseline  as  a weighted  average  of  the  P-R  and  T-P  intervals, 
and  to  make  the  baseline  potential  be  0 mV,  which  will  simplify  sub- 
sequent amplitude  measurements. 

The  traditional  technique  for  straightening  and  shifting  the 
baseline  to  zero  has  been  to  use  a high-pass  filter,  which  removes  the 
low-frequency  baseline  (<0.2  Hz)  while  passing  the  heart  signal  (0.2  Hz 
< heart  signal  < 100  Hz)  [ 7].  This  approach  has  worked  fairly  well  in 
the  past  when  the  subsequent  analysis  was  performed  by  physicians; 
however,  filters  are  not  perfect,  and  the  resulting  distortion  of  the 
heart  signal  itself  is  easily  detectable  when  either  the  cardiogram 
display  is  enlarged  or  computer  analysis  is  performed.  Because  of  this 
problem  and  because  of  the  availability  of  a digital  computer,  a more 
sophisticated  approach  is  used  that  straightens  and  removes  the  baseline 
without  changing  the  heart  signal  itself. 

The  baseline -removal  strategy  is  based  on  the  idea  of  removing  the 
P,  T,  and  QRS  waves  from  the  signal  and  approximating  the  baseline  as  a 
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straight  line  fitted  to  the  remaining  data.  The  baseline  estimate  is 
then  subtracted  from  the  original  signal,  as  shown  in  Figure  1.  The 
P,  T,  and  QRS  waves  are  eliminated  from  the  signal  as  follows.  The 
signal  is  partitioned  into  80-msec  segments  (40  samples  at  500  samples/sec), 
and  the  instantaneous  slope  is  checked  at  each  point.  If  at  any  point 
in  an  80-msec  segment  the  instantaneous  slope  is  greater  than  10  mV/sec, 
then  that  segment  of  data  is  eliminated.  If  it  is  not,  the  value  of  the 
signal  at  the  end  of  the  80-msec  segment  is  recorded  and  used  to  cal- 
culate the  baseline.  For  time  segments  which  contain  either  QRS  com- 
plexes or  T and  P waves,  the  value  of  the  signal  saved  for  the  previous 
time  segment  is  recorded  and  used  in  estimating  the  baseline. 

After  11  values  of  the  signal  are  recorded  (800  msec) , the  baseline 
is  estimated  at  the  midpoint  of  the  800-msec  interval  of  data  by  fitting 
a straight  line  to  the  11  recorded  points.  The  oldest  data  point  is 
then  discarded,  a new  80-msec  segment  of  data  is  checked,  and  either  a 
new  poi.nt  is  recorded  or  the  last  point  is  shifted  down.  The  estimate 
of  the  baseline  at  the  midpoint  of  the  shifted  segment  of  data  is  then 
calculated  and  the  process  repeated.  For  the  80  msec  between  the 
estimates,  linear  interpolation  is  used  to  estimate  the  baseline. 


Figure  1.  baseline  removal. 
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The  straight  line  that  is  fitted  to  the  11  points  (80-msec  segments, 
800-msec  total)  is  of  the  form 

x(k)  = b + ks 

where  k refers  to  the  11  samples  (k  = -5,  -4,  -3,  -2,  -1,  0,  1,  2,  3,  4,  5) 

0 80  160  240  320  400  480  560  640  720  800  Time 

(msec) 


-5-4  -3  -2  -10  1 2 3 4 5 k 

As  indicated  in  the  Slope  Estimation  section,  this  particular 
indexing  of  time  makes  the  estimate  of  the  baseline  at  the  midpoint  of 
the  800-msec  window  of  data  be 
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b = ]Tx(k)  (5) 

k =-5 

To  better  understand  how  the  baseline-removal  logic  is  implemented, 
consider  thi  following  typical  case. 

Situation  1 

a.  Array  W2  contains  data  with  baseline  removed  and  is  600 
samples  long. 


b.  Array  ST  contains  the  latest  data  received  and  is  200  samples 


long. 


c.  Array  SAV  contains  the  last  40  samples  read  into  W2. 


"599! 


600  ST  799 


560  SAV  599  j 

In  the  above  diagram,  W2  is  an  array  of  600  data  points 


fc0  fcl  t2  t3 


W2  (0) 
W2  (1) 
W2  (2) 


x (tQ) 

X (t  ) 
X (t^) 


W2  (598) 
W2  (599) 


x !>98! 
x "599’ 


ST  is  an  array  of  200  data  points 


ST  (0) 
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and  SW  is  an  array  of  40  data  points 


SAV  (0) 
SAV  (1) 


x ‘W 


SAV  (38) 
SAV  (39) 
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X (t599) 


Step  1 

Move  SAV  ahead  40  samples: 


Move  ST  ahead  40  samples : 


I 600  ST  799  | 
I 600  SAV  639  I 


W2  599  600  SAV  639 


ST  8 


Step  3 

Using  data  through  839th  sample,  calculate  the  baseline  at  the  639th 


1 


Step  4 

Using  present  and  previous  (at  599th  sample)  baseline  estimates  and 
linear  interpolation,  remove  the  baseline  from  data  in  SAV  and  move  W2 
ahead  40  samples 

140  W2  639~1  | 640  ST  839  | 

| 600  SAV  639  | 


The  11  data  points  needed  to  estimate  the  baseline  at  the  639th 
sample  are  stored  in  the  arrays  WYR,  ZYREM,  and  SAV.  ZYREM  is  an  array 
of  6 data  points,  and  WYR  an  array  of  5.  The  latter  two  arrays  are 
operated  as  a pushdown  register.  At  the  beginning  of  the  last  example, 
every  40th  point,  beginning  with  the  399th  sample,  was  stored  in  WYR  and 


ZYREM. 


WYR (5) 
WYR  (4 ) 
WYR(3) 
WYR (2) 
WYR (1) 
ZYREM (6) 
ZYREM (5) 
ZYREM (4) 
ZYREM (3) 
ZYREM (2) 
ZYREM (1) 


x 

x 

x 

x 

X 

X 

X 

X 

X 

X 

X 


(t799* 

(t759> 
(t  ) 

t719 

t679 

t639 

t5" 
559 
(t  ) 
519' 
(t  ) 

t479 

t439 

399 


/ 


\ 


(9) 


In  terms  of  the  previous  diagrams,  we  have  the  following  situation. 


Situation  2 

Same  as  situation  1 except  the  arrays  ZYREM  and  WYR  have  been  added. 
Also  the  baseline  estimate  is  known  a*-  the  599th  sample. 


W2 


399 

439 

479 

519 

559 

5991 

T 

ZYREM (1) 

f 

ZYREM (2) 

t 

ZYREM (3) 

1 

ZYREM (4) 

j 

ZYREM (5) 

1 560  5|9 | 

ZYREM (6) 

| 600 

T 

T 

719 

7 

799  I 

1 1 

A 

WYR ( 1 ) 

T 

WYR(2) 

4 

WYR (3) 

r 

WYR (4) 

T 

WYR(5) 

11 
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Step  1 

Move  SAV  ahead  40  samples  (see  situation,  step  1) : 


Step  2 

Store  ZVREM(l)  and  move  ST,  WYR,  and  ZYREM  ahead  40  samples. 


W2 


599] 


399  439  479  519  559  59< 

f f I 1 ( 

2SflVl  ZVrJ.11)  ZVrJ(3)  ZYmHW  ZVrJ.5) 


SAV 


ST 


600 


T 

ZYREM (6) 


|9|  | 640  6^9  719  7^9  799  839| 

MVR  111  ^ WYR  f ^ WYR 


'•WRll)  WVR  (2)  "W31  «i(«)  WY“,5) 


Step  3 

Using  data  through  the  839th  sample,  calculate  the  baseline  at  the  639th 
sample . 

Step  4 

Using  present  and  previous  baseline  estimates  and  linear  interpolation, 
remove  the  baseline  from  data  in  SAV  and  move  W2  ahead  40  samples. 


40 


439 

f 

ZYREM (1) 


479 

* 


ZYREM (2) 


W2 

519 

* 


559  599 

1 1 


ZYREM (3) 


ZYREM (4) 


ZYREM (5) 


600 

639 1 

f 600 

6391 

SAV  j 

ZYREM (6) 

719  759  799  8391 

t f t | t 

ni  1 WYR  (3)  1 WYR  ( ! 


|640  679 

WYR ( 1) 


WYR (2 ) 


WYR (4) 


WYR(5) 


At  the  third  step,  the  data  needed  to  estimate  the  baseline  (see 
equation  5)  at  the  639th  sample  are  stored  in  the  arrays  ZYREM  and  WYR: 


b (639) 


_1 

11 


5 

£ [zYREM(i)  + WYR(i)]  + ZYREM  (6) 
i=l 


(10) 


The  data  needed  to  estimate  the  baseline  at  the  599th  sample  (40  samples 
ago)  were  stored  in  the  arrays  ZYREM  and  WYR  (see  equation  9) 


b (599) 


_1 

11 


5 

£ [ZYREM  x(i)  + WYR_1(i)]  +ZYREM_1(6) 


The  subscript  "-1"  indicates  the  values  of  ZYREM  and  WYR  at  step  1. 

In  going  from  step  1 to  step  2,  the  data  points  in  ZYREM  and  WYR  are 
pushed  down,  ^^gg^  is  stored  as  ZSAVl,  and  ^tgjg)  is  added;  now  ZYREM 

and  WYR  contain  the  following  points  (compare  equations  12  and  9) : 


WYR (5 ) 

* (t839) 

WYR (4) 

X (t799} 

WYR (3) 

x (t '”) 

WYR ( 2 ) 

- :Q 

WYR (1) 

X (t6795 

ZYREM (6) 

X (t639} 

ZYREM (5) 

X (t599> 

ZYREM (4) 

x (t|”l 

ZYREM (3) 

x (t519) 

ZYREM (2) 

X *”479 

ZYREM (1) 

_x  (t439lJ 

The  new  sample,  x(t  ),  is  stored  in  WYR<5)  if,  and  only  if,  the 
o jy 

data  from  t_on  to  t pass  certain  tests  indicating  that  neither  a P, 
799  o39 

T,  nor  QRS  wave  occurs  during  that  80-msec  time  interval.  Finally, 
equations  10  through  12  yield  the  following  recursive  formula  for 
estimating  the  baseline: 

b(k)  = b (k-40)  + yy  (zYREM(l)  - ZSAVl)  (13) 

The  baseline  estimate  is  initialized  via  equation  10,  and  determined  by 
equation  13.  Equation  10  requires  10  additions  and  1 division,  while 
equation  13  requires  only  2 additions  and  1 division. 

There  are  two  tests  for  determining  whether  or  not  an  80-msec 
interval  of  data  contains  a P,  QRS,  or  T wave:  (1)  The  maximum  negative 
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slope  of  the  R wave  in  the  x lead  is  used  as  a fiducial  point  in  per- 
forming signal  averaging  (see  Digital  Filter  section).  Therefore,  since 
the  fiducial  point  is  available,  S-T  intervals  are  eliminated  by  rejecting 
from  the  baseline  estimation  all  data  occurring  during  the  100  msec 
following  a fiducial  point.  (2)  The  instantaneous  slope  is  checked  at 
every  point  in  the  bO-msec  interval.  If  that  slope  exceeds  6 mV/sec  at 
any  point  in  the  interval,  it  is  concluded  that  either  a P,  QRS,  or  T 
wave  occurred  during  that  interval.  The  6-mV/sec  threshold  was  deter- 
mined experimentally  and  is  valid  only  with  the  7-point  slope  estimate 
of  equation  4,  which  can  be  rewritten  as 

Q71 

(14) 

Each  of  the  three  terms  enclosed  by  brackets  in  equation  14  is  an 
estimate  of  the  slope  of  the  waveform.  The  actual  estimate  of  the  slope 
at  the  kth  data  point  is  merely  the  weighted  average  of  the  three 
estimates , with  the  weightings  proportional  to  the  time  intervals 
involved.  For  instance,  the  first  difference,  stretching  over  six 
samples  (12  msec),  is  given  the  heaviest  weighting  (.643).  These 
weightings  coupled  with  the  6-mV/sec  threshold  are  important  in  dis- 
tinguishing the  P,  QRS,  and  T waves  from  signal  noise. 


jx(ktl)-x(k-l)j 


, .,|x(k+3)-x(k-3)1 
s(k)  = .643[ 1 -J 


The  actual  version  of  equation  4 used  in  the  computer  program  is 
recursive  in  nature  and  requires  fewer  additions  and  multiplications. 
The  version  used  follows  from  comparing  equation  15 


s (k  + 1)  = .643 


[ 


x(k+4) -x(k-2) 


J +.286  | 


x (k+3) -x (k-1) 


.071  [ 


115) 


with  equation  14,  which  yields 

s(k+l)  = s (k)  - { £ (k)+3  [x(k+4)+x(k-3)j  } /28  (16) 

where  £ is  the  solution  of 


I (k+1)  = £ (k)  + x (k+4)  -x  (k—  2 ) (17) 

The  slope  estimate  is  initiated  via  equation  14  and  then  is  propagated 
via  equation  16.  The  sum  £ , used  in  equation  16,  is  initiated  via 

£ (k)  = x(k+3)  + x(k+2)  + x(k+l)  + x(k)  + x(k-l)  + x(k-2) 


and  is  then  propagated  via  equation  17. 


Up  to  this  point,  we  have  presented  an  algorithm  which  estimates 
values  of  the  baseline  at  every  40th  sample  (every  HO  msec). 


Time  (msec) 


To  subtract  the  nonzero  baseline  from  the  original  signal,  we  need 
an  estimate  of  the  baseline  at  every  sample  (every  2 msec) . This  is 
generated  by  fitting  a straight  line  between  each  pair  of  points. 


Time  (msec) 


Finally,  the  signal  baseline  is  shifted  to  zero  by  subtracting  the  fine 
baseline  estimate  from  the  original  cardiogram. 


Fiducial  Detection  and  Averaging 


To  reduce  the  quantity  of  data  to  be  subsequently  analyzed  and  to 
remove  the  noncardiac-generated  noise  that  is  uncorrelated  with  itself 
from  one  beat  to  the  next,  the  digital  filter  will  average  successive 
beats  and  produce  a single  representative  beat  for  subsequent  analysis. 
This  approach  has  been  successfully  used  by  others  [l,  2,  4,  8,  10]. 

The  steps  required  are  to  locate  each  beat,  align  them  temporally,  add 
the  beats  together,  and  divide  the  resulting  sum  by  the  number  of  beats 
used.  The  point  used  in  locating  each  beat  is  the  point  on  the  x lead 
corresponding  to  the  maximum  absolute  value  of  its  slope  or  time  deriv- 
ative. This  fiducial  point  has  been  successfully  used  by  others  [l,  4] 
and  is  relatively  insensitive  to  baseline  disturbances. 


The  instantaneous  time  derivative  estimate  used  in  the  baseline- 
removal  algorithm  (equations  16  and  17)  is  also  used  to  detect  fiducial 
points.  The  algorithm  for  detecting  these  points  first  determines  the 
maximum  absolute  value  of  the  slope  of  the  x lead  in  the  first  1.6  sec 
of  data.  This  maximum  value  is  used  to  determine  a threshold,  which  is 
set  at  60%  of  the  max  value.  After  the  threshold  is  initialized,  the 
slope  is  calculated  and  compared  to  the  threshold  as  the  data  are 
received.  After  the  threshold  is  exceeded,  the  point  of  maximum  slope 
is  used  as  the  fiducial  point;  the  threshold  is  then  reset  to  60%  of  the 
maximum  value  just  determined,  and  the  fiducial  point  detector  is 
deactivated  for  the  next  400  msec  to  prevent  double-peaked  R waves  from 
becoming  a problem. 


Bad-Data  Rejection 

Under  certain  conditions,  it  is  impossible  to  extract  the  actual 
cardiograms  from  the  noise.  Although  it  does  not  cover  all  cases,  the 
criterion  used  to  reject  or  accept  data  is  based  on  the  slope  of  a 
pseudobaseline  estimate.  This  estimate  is  determined  as  was  the  actual 
baseline  estimate  with  the  exception  that  points  occurring  during  the  P, 
QRS , and  T waves  are  used.  Any  time  that  the  slope  of  this  straight 
line  exceeds  2 mV/sec,  the  data  are  not  used  in  calculating  the  average 
waveform. 

The  pushdown  register  used  to  store  the  data  needed  to  calculate 
the  pseudobaseline  (11  successive  points,  each  80  msec  apart)  is  com- 
posed of  the  arrays  YREM  and  ST. 


YREM(l) 
YREM (2) 


YREM (6) 
ST (39) 
ST (79) 


[ST( 199)1 

With  N = 11,  equation  4 can  be  manipulated  (as  in  the  Baseline-Removal 
section)  to  obtain 

s(k-t-l)  = s(k)  - { I(k)+5  [ x (k+6)  + x(k-5)]  } /110  (18) 

where 

I(k+1)  = I (k)  + x (k+6)  - (k-4)  (19) 
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Computer  Program  Flowchart 

A list  of  the  key  variables  and  a flowchart  of  the  digital  filter 

program  follow. 

Variables 

TREM  - Variable  used  with  the  peak-detection  scheme. 

LC  - Counter  activated  after  fiducial  has  been  located.  Used  to  disable 
fiducial  detection  for  65  samples. 

ITR  - Flag  which  is  set  to  1 when  instantaneous  slope  is  less  than 
threshold  (fiducial  detector  is  enabled) . Reset  to  0 when 
fiducial  is  detected. 

K - Dummy  variable  used  in  initialization  of  average  slope  estimator. 

JF  - Counter  activated  when  average  slope  exceeds  2 mV/sec.  Stops 
averaging  for  100  samples. 

IR(10,3)  - Counters  activated  when  fiducial  point  is  detected. 

IRC  - Counts  the  number  of  fiducials  detected. 

ZYSUM ( 3 ) - Sum  of  11  points,  once  every  20  samples  and  with  QRS  points 
excluded. 

YSUM(3)  - Sum  of  11  points,  once  every  20  samples  with  no  points  excluded. 

TYSUM<3)  - Numerator  of  least-squares  estimate  of  slope,  extends  over 
200  samples. 

NW(3)  - Counter  used  to  record  number  of  beats  used  to  form  average 
beat. 

NR<3)  - Counters  used  to  determine  number  of  samples  between  average 
waveform  update. 

JYF  ( 3 ) - Flags  set  when  instantaneous  slope  exceeds  10  mV/sec.  Used  to 
exclude  data  from  baseline  estimate. 

W2  (299,3)  - Contains  waveforms  with  baseline  removed. 

AV(299,3)  - Contains  averaged  waveforms. 

THRESH  - Threshold  used  in  peak-detection  logic. 

TXS(3)  - Instantaneous  V slope  estimate  for  each  lead. 
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:’.YREM(6,3)  - Array  used  with  baseline  estimator.  Points  with  bad 
baselines  are  removed. 

IRC  - Number  of  beats  detected. 

YREM (6,3)  - Array  used  with  average-slope  estimator.  Contains  all 
points — none  are  rejected. 


/ 
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INITIALIZE  DIGITAL  FILTER 


INITIALIZE  JRS 


INITIALIZE 
INSTANTANEOUS 
SLOPE 
ESTIMATOR 
(TXS,  DSUM) 


ENTER  NEW  DATA 


TXS(LD)  = 4-JRS(l  ,LD)  + 3‘M(LD)  + TXS(LD)  - DSUM(LD) 
DSUM<LD)  = DSUM(LD)  - JRS(l.LD)  + M(LD) 

r 


CALCULATE 

INSTANTANEOUS 

SLOPE 


JRS<7,LO»=M(LD)| 


INCREMENT 
v INSTANTANEOUS 
SLOPE  SHIFT 
REGISTER 


CHECK  TO  SEE 
IF  THIS  80  msec 
CORRESPONDS 
TO  A FLAT  PORTION 


OF  THE  CURVE 


YREM(I,LD)=M(LD) 


ZYREM(I.LD)=ZYREM{M.LD) 

JYF(LD>=0 


ZYREM(1.LD)=M(LD) 


ZYSUM(LD)=ZYSUM(LD)+ZYREM(I.LD) 
YSUM<LD^YSUM(LD)+YREM(I,LD) 
TYSUM(LDhTYSUM(LD)-KFLOAT  (K-240))  M(LD) 


i 


WYR(*.LD)=ZYREM(6.LD) 


1 1 NO 

LD=LD+1  f-*— 


| INITIALIZE  BASE 
I LINE  AND  LONG 
rTERM  SLOPE 
ESTIMATORS 


INITIALIZE  BASELINE 
AND  LONG-TERM 
SLOPE  ESTIMATES 


INITIALIZE 
COUNTERS 
AND  TIME 


STORE  OLDEST 
POINTS  USED  IN 
BASELINE  AND 
SLOPE  ESTIMATES 


STORE  OLDEST 
42  POINTS  IN 
ST  IN  SAV 
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t 


1 


ZYREM(1.LD)=ZYREM  (I+l.LD) 
YREM(J.LD)=YREM(I+1,LD) 
ZYREM(6.LD^WYR(1.LD) 
YREM(6.1D>=$T(JST,LD) 


SHIFT  BASELINE 
AND  SLOPE  - 
ESTIMATE 
COUNTERS 


START  OF  BLOCK  TO  READ 
IN  8 1 msec 
WORTH  OF  DATA 


TXS(LD)=4-JRS(l.LDH-3-V{LDKTXS(LD)-DSUM(LD) 
DSUM(LD)=DSUM(LD)-JRS(1.LD}+M(LD) 
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IDF18 


L 
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[IDF  191 


IDF  131 


, NO  / \ 

I=Wk-/l=39, 


lST=MOD(IST+40.200) 


END  OF  80-msec  BLOCK 


IDF20) 


WYR(1.LD)=WYR(I+1,LD) 


NO  /\ 

1=1+1  h < 1=4 


JYF(LD)=1 


JYF(LD)=0 


INCREMENT  BASELINE 
ESTIMATOR  STORAGE 


WYR(5.LD)=M(LD) 


TYSUM(LD^TYSUM(LDH-2  4E2  SAV1(LD) 
+2E2’M(LD)-4ErYSUM(LD) 
ZYSUM(LD)=ZYSUM(LD)-ZSAV1(LDH-WYR(5.LD) 
YSUM(LD)=YSUM(LO)-SAVl(LD}fM(LD) 
B(LD)=ZYSUM(LD)/1 1 
SDJF(LDHB(LD)-Bl(LD)y40 
S(LD)=TYSUM(LD) 


CALCULATE 
BASELINE  AND 
SLOPE  ESTIMATE 

-*-/lDF2l\ 
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CLEANUP  BASELINE 
FOR  DATA 
STORED  IN  SAV 


1 W2=MOD<I  W2  +40.600) 


« 

STORE  BASELINE 
ESTIMATE 
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Program  Verification 


The  program  logic  was  verified  both  at  USAFSAM  and  at  Draper 
Laboratory.  The  verification  tests  can  be  grouped  as 

(1)  Baseline-removal  test 

(2)  Piducial-point-detection  test 

(3)  Averaging  test 

(4)  Data-good,  Data-bad 

At  USAFSAM,  the  tests  were  performed  as  follows.  A nurses'  training  aid 
was  used  to  generate  the  reference  heart  signal  represented  in  Figure  3. 


Figure  3.  The  reference  cardiogram. 


The  reference  signal  was  then  fed  into  the  system.  Results  are 
shown  in  Figure  4:  The  top  trace  (a)  of  Figure  4 corresponds  to  the 

original  reference  cardiogram;  the  middle  trace  (b)  is  the  output  of  the 
baseline-removal  algorithm;  and  the  bottom  trace  (c)  is  the  difference 
between  the  output  of  the  algorithm  and  its  input.  The  slight  d.c.  bias 
has  been  removed  so  that  no  unwanted  signals  have  been  added.  A 3-in 
(7.62-cm)  vertical  line  is  printed  to  signify  where  fiducial  points  have 
been  detected.  The  averaging  logic  writes  a 0.31-in  (0.79-cm)  vertical 
line  on  the  middle  trace,  immediately  after  a complete  beat  has  been 
located  and  added  to  the  average.  In  this  test,  all  17  beats  were 
detected  and  used  to  form  the  averaged  beat  that  is  displayed  in  Figure  5 

Next,  the  reference  cardiogram  was  corrupted  by  adding  an  additional 
sine  wave  baseline.  Three  different  frequencies  were  used  (0.079  Hz, 
0.158  Hz,  and  0.5  Hz),  and  for  each  frequency  three  different  amplitudes 
were  used  (the  lowest,  not  to  trip  the  data-good  data-bad  logic;  the 
middle,  to  trip  it  occasionally;  and  the  highest,  to  trip  it  every 
time . ) 

The  results  for  the  0.5-Hz  sine  wave  are  given  first.  Figures  6 
and  9(a)  show  the  results  of  choosing  a 0.5656-mV  rms  amplitude  to  make 
the  2-mV/sec  data-bad  limit  always  be  exceeded;  no  beats  were  averaged. 
All  fiducial  points  were  detected  and  labeled,  and  the  0.5-Hz  sine  wave 
was  attenuated  =»20%. 


Figure  5.  The  average  reference  cardiogram. 


Figures  7 and  9(b)  show  the  system's  performance  with  the  same 
frequency  noise  (0.5  Hz),  but  at  an  amplitude  (0.2828-mV  rms)  carefully 
chosen  to  make  the  data-bad  limit  be  exceeded  at  times  and  not  exceeded 
at  others.  As  expected,  some  beats  were  used  in  calculating  the  average 
beat,  while  others  were  not;  also,  the  fiducial  point  was  detected  for 
every  beat.  Finally,  although  the  0.5-Hz  contamination  was  slightly 
attenuated  (possibly  20%),  the  averaged  beat  (Fig.  9(b))  was  slightly 
tilted. 

To  finish  the  0.5-Hz  test,  an  amplitude  was  chosen  that  would  not 
exceed  the  2-mV/sec  data-bad  logic.  The  results  for  this  amplitude, 
0.1414-mV  rms,  are  shown  in  Figures  8 and  9(c).  All  fiducials  were 
detected,  and  as  expected,  all  beats  were  used  to  calculate  the  averaged 
beat.  The  contaminating  sine  wave  was  attenuated  20%  by  the  filter,  and 
because  of  the  lower  amplitude  of  the  sine  wave,  the  averaged  beat  was 
less  tilted. 

Next,  the  frequency  of  the  noise  signal  was  decreased  to  0.158  Hz, 
and  the  same  sequence  of  amplitudes  (high,  to  always  trigger  data-bad 
test;  medium,  to  trigger  it  some  time  but  not  always;  low,  to  never 
trigger  data-bad  test)  were  used  to  generate  baseline  noises.  The 
results  obtained  with  the  high- amplitude  (2.121-mV  rms)  noise  are  given 
in  Figures  10  and  13(a).  All  fiducials  were  identified,  and  as  expected, 
the  data-bad  criterion  was  always  exceeded;  the  average  waveform  was  not 
calculated.  The  0.158-Hz  contamination  was  attenuated  80%.  Then  the 
amplitude  of  the  sine  wave  contamination  was  decreased  to  1.061-mV 
rms,  a value  that  should  make  the  data-bad  limit  be  exceeded  at  times 
and  not  at  others;  this  was  the  case,  as  shown  in  Figures  11  and  13(b). 
Finally,  the  amplitude  of  the  sine  wave  contamination  was  decreased  to 
0.55146-mV  rms,  a level  that  should  make  all  the  data  acceptable.  As 
indicated  by  Figures  12  and  13(c),  all  fiducials  were  detected  and  all 
beats  were  used  in  determining  the  averaged  beat. 

The  same  procedure  was  repeated  with  a 0.079-Hz  sine  wave  (a  more 
realistic  baseline).  The  results  (using  4.242-,  2.121-,  and  1.06-mV 
rms) , which  were  the  same  as  before,  are  shown  in  Figures  14  through  17. 


Figure  7.  Test  results  for  reference  signal  corrupted  by  a 0.5-Hz,  0.2828-n\V  rms  sine  wave. 
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Figure  10.  Test  results  for  reference  signal  corrupted  by  a 0.158-Hz,  2.121-mV  rms  sine  wave. 


Figure  ’l.  Test  results  for  reference  signal  corrupted  by  a 0.158-Hz,  1.061-mV  rms  sine  wave. 
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Figure  14.  Test  results  for  reference  signal  corrupted  by  a 0.079-Hz,  4.242-mV  rms  sine  wave 


Figure  15.  Test  results  for  reference  signal  corrupted  by  a 0.079-Hz,  2.121-mV  nns  sine  wave. 


Several  tests  were  made  at  the  Charles  Stark  Draper  Laboratory  to 
precisely  examine  the  response  of  the  baseline-removal  algorithm  to 
idealized  high-energy  waveforms.  The  objectives  were  to  determine  if 
any  computer-generated  artifacts  occurred  and  to  determine  the  opera- 
tional limits  of  the  algorithm.  Three  types  of  waveforms  were  used  in 
the  tests:  (1)  triangular  waves,  (2)  square  waves,  and  (3)  pure  steps. 

Results  obtained  using  a very-high-energy  triangle  waveform  are 
seen  in  Figure  18.  The  triangles  have  a 4-mV  amplitude  and  a sinusoidal 
baseline  with  a 2-mV  amplitude  and  0.16-Hz  frequency.  With  these  values 
the  maximum  baseline  slope  is  2.01  mV/sec,  just  exceeding  the  2-mV/sec 
threshold  for  bad-data  rejection;  however,  no  bad-data  tests  were  used 
here — the  objective  was  to  evaluate  the  baseline-removal  algorithm.  The 
underlying  heart  rate  is  100  beats/sec,  which  is  borderline  tachycardia. 
Higher  heart  rates  were  used  to  insure  that  QRS  complexes  in  close 
proximity  would  not  lead  to  problems  in  reconstruction.  Figure  18 
demonstrates  good  waveform  reconstruction.  The  only  computer-generated 
artifact  is  a small  residual  baseline  which  tends  to  average  out  in 
time.  Slope  estimate  is  affected  by  the  large  triangular  waveforms,  but 
is  centered  around  the  correct  value  and  exhibits  no  phase  shift. 


Next,  a series  of  square  waves  were  made  using  amplitudes  of  up  to 
4 mV  and  durations  of  up  to  160  msec.  In  Figure  19,  the  results  obtained 
using  a 4-mV  square  wave  of  80-msec  duration,  superimposed  on  the  same 
sinusoidal  baseline,  are  shown.  The  baseline  removal  is  excellent,  with 
only  a small  residual  artifact  that  depends  upon  the  baseline  remaining. 
The  slope  estimate,  while  affected  by  the  square  waves,  is  correctly 
centered  and  exhibits  zero  phase  shift.  Square  waves  of  duration 
greater  than  80  msec  could  not  be  handled.  If  duration  is  greater,  two 
samples  may  occur  within  the  square  wave,  indicating  zero  slope. 

Although  provision  could  easily  be  made  for  handling  wider  square  waves, 
this  does  not  appear  necessary  since  waves  are  nonphysiological  in 
origin. 


Figura  19.  Baseline-removal  evaluation,  square  QRS  wave. 


Finally,  tests  were  made  to  determine  the  response  to  sudden  step 
changes  in  the  baseline.  Figure  20  shows  the  results  of  a typical  test. 
First  a jump  of  0 to  1 mV  was  used,  with  a subsequent  jump  of  -2  to 
-1  mV.  The  corresponding  baseline  and  slope  estimates  are  shown  in  the 
figure.  The  effects  of  the  step  on  the  estimates  are  spread  out  evenly 
on  both  sides  of  the  jump,  owing  to  the  symmetry  of  the  moving  window. 
There  is  a baseline  error  within  220  samples  (440  msec)  of  the  jump,  and 
a slope  error  within  200  samples  (400  msec)  of  the  jump.  The  maximum 
baseline  error  is  0.5  mV/sec  per  mV,  and  the  maximum  slope  error  is 
about  2 mV/sec  per  mV.  The  results  of  using  the  idealized  waveforms  of 
Figures  18-20  indicate  that  the  algorithm  operates  satisfactorily  in 
these  cases  and  introduces  only  small  artifacts  that  tend  to  average  out 
with  time. 
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Figure  20.  Baseline-removal  evaluation,  step  baseline  changes 


OUTPUT  DATA  FORMAT!? 

Three  types  of  data  output  are  generated  by  the  processing  system. 
All  files  are  on  0.5-in  (1.77-cm)  9-track  magnetic  tape  recorded  at  800 
BPI.  Two  Of  these  are  the  standard  RDOS  formatted  files,  and  the  third 
is  written  for  input  to  an  IBM  360/65  for  further  analysis  of  the 
vectors  ard  amplitudes  which  cardiologists  use  to  identify  coronary 
problems  [ 1 1 ] . Dogs  of  all  tapes  are  created  on  the  line  printer  as  the 
data  are  stored  on  tape. 


- 


RDOS  File  Formats 


Both  of  the  RDOS  files  use  the  standard  system  of  512  bytes  of  data 
and  2 bytes  of  file  identification.  These  files  are  buffered  and  may  be 
referenced  an  element  at  a time.  The  data  tape  produced  by  the  data 
acquisition  or  digitizing  program  consists  of  a 6-byte  identifier, 
normally  the  vectorcardiogram  number,  right- justified  with  leading 
zeros.  This  is  followed  by  three  16-bit  integers  containing  the  age, 
height,  and  weight  of  the  subject.  Then  follows  2250  16-bit  integers  of 
calibration  pulses,  750  per  channel,  written  X^,  Y , Z1 , , 1^, 

X750'  Y750'  Z750'  FlnaHy  18,000  16-bit  integers  of  vector- 
cardiogram data  are  written  in  the  same  manner  as  the  calibration  data. 
This  is  followed  by  an  end-of-file  mark.  A segment  of  data  is  checked 

for  overscaling  before  it  is  transferred  to  magnetic  tape.  If  over- 
scaling of  calibration  data  has  occurred,  a message  is  written  on  the 

line  printer  and  the  control  terminal  and  a restart  occurs;  if  not,  the 

data  is  written  on  tape.  If  overscaling  occurs  on  the  vectorcardiogram, 
the  message  is  written  again  and  the  operator  can  restart. 


The  output  of  the  processing  program  described  in  the  Digital 
Filter  section  contains  the  same  identification,  age,  height,  and 
weight  information  as  the  digitizing  program.  This  is  followed  by  three 
real  values  which  are  the  scale  factors  for  X,  Y,  and  Z leads  respec- 
tively. To  transform  the  filtered  vectorcardiogram  data  to  millivolts, 
each  value  is  divided  by  the  scale  factor  for  that  lead.  After  the 
scale  factors,  come  16,800  16-bi.t  integers  of  the  filtered  vectorcardio- 
gram in  Xr  Yr  Zx,  X2,  Y2,  Z2,  . . . X560Q,  Y56qo,  Z^  sequence. 

These  data  are  followed  by  a sentinel  integer  of  32767  for  each  lead. 
Then  comes  one  real  value  which  is  the  heart  rate  for  this  segment. 

Next,  1800  real  values,  representing  the  averaged  and  calibrated  wave- 
forms, are  written  in  X,,  Y, , Z, , X„,  Y„,  Z^,  . . . X„„„ , Y, , Z, 


sequence . 


1’  1'  1'  2'  2'  2'  600'  600'  “600 
The  patient  record  is  closed  with  an  end-of-file  mark. 


IBM  360/65  Tape  Formats  1 

The  data  for  input  to  the  IBM  360/65  are  the  same  as  the  output  of 
the  processing  program  but  have  been  written  in  six  blocks  of  7200 
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bytes  each.  The  first  block  contains  the  identification  information  in 
ASCII  characters:  the  age,  height,  weight,  scale  factors,  and  3588 

vectorcardiogram  values.  Blocks  2,  3,  4,  and  5 contain  the  remainder  of 
the  -'ectorcardiogram,  the  sentinel,  and  the  heart  rate.  Block  6 con- 
tains the  1800  real  values  of  the  averaged  waveform. 


DIGITIZING-SYSTEM  OPERATION 

To  operate  the  programs  used  to  acquire  and  process  cardiograms, 
the  operator  must  be  familiar  with  the  Data  General  RDOS . This  system 
is  described  in  Data  General's  manual  entitled  "RDOS,  Real-Time  Disk 
Operating  System." 

Once  familiar  with  this  system,  the  operator  can  use  the  programs 
to  obtain  data  and  request  actions  through  a series  of  questions  and 
messages  on  the  operator's  console.  Certain  pieces  of  hardware  must  be 
operational  and  identified  to  the  system  for  the  programs  to  operate 
properly. 


Required  Programs 

Three  programs  are  required  to  complete  the  process  from  data 
acquisition  to  IBM- 360-compatible  tape  generation.  The  first  program, 
DRADIG,  acquires  the  data  and  creates  a magnetic  tape  file  for  each 
subject.  The  second  program,  F5DRA,  processes  the  output  of  the  first 
program  by  determining  calibration  factors,  filtering  the  vectorcardio- 
gram data,  and  producing,  for  each  lead,  an  average  beat  calibrated  to 
millivolts.  This  program  also  produces  an  output  file  for  each  program. 
The  third  program,  T0360,  converts  the  output  of  F5DRA  to  IBM-360- 
compatible  tape. 

Required  Hardware 

The  basic  hardware  requirements  for  the  current  system  are  two 
9- track  magnetic  tapes  (800  BPI) , a fixed-head  NOVADISK  for  temporary 
data  storage  while  digitizing,  a graphics  terminal  connected  to  channel  1 
of  the  multiplexer,  and  at  least  three  channels  of  analog-to-digital 
(A/D)  converters.  Also  required  are  a 32K  NOVA  minicomputer  with 
floating-point  processor  and  hardware  multiply  and  divide,  and  an  analog 
magnetic  tape  transport  with  the  X channel  connected  to  A/D  channel  0,  Y 
to  channel  1,  and  Z to  channel  2. 


Operating  Procedures 

DRADIG — The  data  acquisition  program,  DRADIG,  obtains  all  operator- 
provided  information  and  control  instructions  through  the  operator's 
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console.  A magnetic  tape  must  be  mounted  on  unit  zero  and  loaded.  The 
tape  is  identified  to  the  system  by  the  operator  typing  INIT  MTO  and 
RETURN.  The  program  can  then  be  started  by  typing  DRADIG  and  RETURN. 

The  program  then  asks  NEW  TAPE,  YES  OR  NO,  and  the  operator  types  the 
proper  response  and  a return.  If  he  types  NO,  the  program  assumes  that 
the  tape  contains  some  data  and  asks  NUMBER  OF  FILES  USED.  The  operator 
responds  with  the  proper  file  number  and  a return.  After  positioning 
the  tape  to  the  proper  file,  the  program  asks  CASE  NUMBER.  The  operator 
responds  with  six  characters,  using  leading  zeros,  and  a return.  The 
program  then  asks  for  age,  height,  and  weight.  Response  to  each  must  be 
an  integer  followed  by  a return. 

The  program  message  DIGITIZING  CALIBRATION — TYPE  RETURN  TO  START 
appears;  and,  when  the  analog  magnetic  tape  is  positioned  at  the  cali- 
bration data  and  running  at  the  proper  speed,  the  operator  types  RETURN. 
Six  seconds  of  calibration  data  are  digitized  by  the  program,  examined 
for  overscale,  and,  if  no  overscaling  is  detected,  spooled  to  magnetic 
tape.  If  overscale  is  detected,  the  bell  on  the  console  rings  and 
numbers  representing  the  number  of  samples  out  of  range  appear  on  the 
console  and  line  printer.  The  message  PAUSE  then  appears,  and  the 
program  waits  for  a return.  After  the  operator  types  RETURN,  the  pro- 
gram asks  for  another  case  number  and  continues  as  before. 

The  program  message  DIGITIZING  DATA — TYPE  RETURN  TO  START  appears; 
and , when  the  analog  tape  is  positioned  to  the  data  segment  required  and 
running  at  the  proper  speed,  the  operator  types  RETURN.  These  data  are 
also  examined  for  overscale,  and  transferred  to  tape  if  none  occurs.  If 
overscale  occurs,  the  same  procedure  applies  as  for  calibration  over- 
scale. After  the  data  have  been  transferred  to  tape,  the  program  asks 
for  another  case  number.  If  the  operator  wishes  to  terminate  the  run, 
he  types  DONE  and  RETURN.  If  he  wishes  to  continue,  he  supplies  a new 
case  number  and  a return. 

F5DRA — The  digital  filtering  and  averaging  program  obtains  all 
control  information  from  the  operator's  console  and  does  any  graphic 
output  on  the  graphic  terminal  connected  to  multiplexer  channel  1. 

Under  normal  automatic  processing,  a data  tape  must  be  mounted  on  one 
tape  transport  and  loaded,  and  an  output  tape  mounted  on  a second 
transport  and  loaded.  Both  must  be  identified  to  the  system  with  the 
INIT  instruction. 

To  load  the  program,  the  operator  types  F5DRA  and  RETURN.  When  the 
program  message  PLOT,  O-NONE,  1-X,  2-Y,  3-Z  appears,  the  operator  types 
the  proper  value  to  eliminate  plotting  or  selects  the  lead  to  be 
plotted,  and  a return.  The  program  message  MODE;  AUTOMATIC  OR  SINGLE 
appears;  and  the  operator  types  either  AUTO  or  SINGLE,  and  RETURN.  A 
response  of  SINGLE  will  cause  the  program  to  terminate  after  processing 
only  one  file;  AUTO  will  cause  the  program  to  process  sequential  files 
from  the  input  magnetic  tape  until  end  of  tape  is  encountered.  If  files 
from  the  disk  are  to  be  processed,  the  SINGLE  mode  is  selected;  and  the 
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program  then  asks  FIRST  INPUT  FILE.  If  AUTO  mode  is  selected,  the 
response  is  of  the  form  "MTN:mm" — where  N is  the  unit  number  of  the 
input  tape  transport,  and  mm  is  the  number  of  the  first  file  for  input. 
Next  the  program  asks  FIRST  OUTPUT  FILE.  For  AUTO  mode,  the  same  type 
of  response  is  given  as  was  given  for  the  input  file  in  AUTO.  For 
SINGLE  mode,  each  response  is  the  6-character  file  name. 

The  program  then  proceeds  to  filter,  scale,  and  average  the  data, 
using  the  graphic  terminal  to  plot  the  selected  lead  and  the  averaged 
waves.  If  the  hard-copy  unit  is  on,  it  also  makes  hard  copy  of  each 
screen  full  of  data.  The  filtered  and  averaged  data  are  written  to  the 
output  file.  Individual  R-R  intervals,  along  with  ratios  of  successive 
intervals,  are  recorded  on  the  line  printer.  The  ratios  are  used  to 
detect  rhythm  disturbances.  The  averaged  waveforms  are  always  plotted. 
The  running  time  of  the  program  is  nearly  doubled  if  plotting  of  a lead 
is  requested. 

T036Q — The  program  to  generate  IBM- 360-compatible  tape  requires  two 
tapes  for  operation.  Both  must  be  loaded  and  identified  to  the  oper- 
ating system  with  the  INIT  instruction.  To  load  the  program,  the 
operator  types  TO360  and  RETURN.  The  queries  and  responses  are  the  same 
as  those  of  the  processing  program  for  mode,  output  and  input  files. 

The  program  sequentially  processes  an  input  tape  and  creates  multiple 
files  on  the  output  tape  for  further  data  processing  on  an  IBM  360/65  at 
San  Antonio  Data  Services  Center. 


CONCLUSION 

A system  for  digitizing  and  filtering  ECG/VCGs  has  been  designed 
and  is  operational  at  USAFSAM.  Over  900  VCGs  have  been  digitized  and 
are  being  used  in  ongoing  research  work,  both  at  USAFSAM  and  under 
contract. 
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